Structural phases of strained LaAlC>3 driven by octahedral tilt instabilities 
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We investigate the effect of epitaxial strain on [001] -oriented LaAlOs using first-principles density 
functional calculations. We find a series of structural phase transitions between states characterized 
by distinct patterns of tilting of the A106 octahedra. By tuning the biaxial strain from compressive 
to tensile, we induce an evolution in the crystal structure in which the tilt axis changes from out- 
of-plane to in-plane, corresponding to space groups H/mcm and Imma. We also study the effect 
of uniaxial relaxation of the usual biaxial epitaxial constraint and explore this as a mechanism for 
selectively stabilizing different patterns of octahedral tilts. 



Oxides in the ABO% perovskite family present a mul- 
titude of functional properties and are widely renowned 
for their potential in technological applications. Con- 
struction of heteroepitaxial thin films is being actively 
explored as a route to further enhance and expand on 
the existing oxide functionalities. The presence of an 
interface between an oxide and a substrate can dramat- 
ically affect material properties, particularly if a film is 
grown coherently so that its in-plane lattice constant is 
forced to match that of the substrate. The resulting het- 
eroepitaxial strain has been credited, for example, with 
dramatically enhancing the ferroelectric polarization and 
Curie temperature in thin film BaTiOs 1 and inducing 
ferroelectricity in usually non-polar SrTiOs. 2 In addition 
to ferroelectric distortions, we increasingly find repercus- 
sions of strain in the patterns of rigid rotations, or tilts, of 
the corner-sharing BOq octahedral units. These changes 
in tilts have been invoked to account for such phenom- 
ena as the strain-dependence of magnetic properties, as 
a result of the corresponding changes in the electronic 
bandwidth*^ In this work, we use density functional the- 
ory to examine the effect of epitaxial strain on the tilt 
instabilities in LaA103 (LAO). LAO is an ideal model 
system in which to isolate the interplay between strain 
and tilting because it exhibits only rotational distortions 
in its ground state with no indications of ferroelectric or 
Jahn- Teller instabilities. 

We are also motivated by a need to understand the 
influence of octahedral rotations in oxide heterostruc- 
tures such as LAO/SrTi03 (STO), which has been re- 
ported to form a highly conductive electron gas at the 
interface in spite of the insulating nature of the two par- 
ent compounds. 5 The propagation of octahedral rotations 
across the interface is likely to be relevant to its electronic 
properties, but microscopic techniques have limited ac- 
cess to the positions of oxygen ions. 6 An improved under- 
standing of the effect of strain on the parent compound 
LAO is therefore highly desirable. 

Here we investigate the effect of biaxial strain in LAO. 
We find that the bulk R3c structure is destabilized by 
the constraint of coherent epitaxy and that remarkably 
small compressive or tensile biaxial strains stabilize pre- 
viously unidentified phases of LAO. In addition, we find 
that uniaxial relaxation of the biaxial epitaxial constraint 
stabilizes a third phase. The primary structural differ- 
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FIG. 1: Energy per 5-atom formula unit of LAO as a function 
of in-plane lattice constant a (constant volume). Energy is 
given relative to unconstrained R3c. Right: Schematics of 
pseudocubic unit cells showing axes of octahedral rotations 
in the three phases. (Color online.) 



ences between these strain-stabilized phases and the par- 
ent phase are found in the patterns of octahedral rota- 
tions, and our analysis of the transition between them 
provides insight into the coupling between strain and ro- 
tations in LAO and similar ABO3 perovskites. 

Calculation details and notation. — For all calculations 
we use the local density approximation (LDA) of den- 
sity functional theory as implemented in the Vienna ab- 
initio simulation package, VASP^ We use the projector 
augmented wave method 8,9 with the default VASP po- 
tentials (La, Al, O) and a plane wave energy cut-off of 
800 eV. We use a 10-atom rhombohedral unit cell with 
a 5x5x5 Appoint sampling. To determine ground state 
structures we relax ionic positions to a force tolerance 
of 1 meV/A. Space group determinations are performed 
with the FINDSYM symmetry analysis software. 10 

Strain is calculated as (a — ao)/ao, where ao is the 
LDA equilibrium lattice parameter. We use the nota- 
tion a, 6, c to denote the pseudocubic (pc) lattice con- 
stants, which correspond to the crystallographic direc- 
tions [100]pc, [010]pc, and [001] pc . Lattice distortions 
comprised of rigid octahedral rotations are described 
by the angles of rotation around the pseudocubic axes, 
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FIG. 2: Transition details as a function of in-plane lattice 
constant a (constant volume). Top panel, energy per 5-atom 
formula unit; middle panel, Al-0 bond distances; lower panel, 
Al-O-Al bond angles. (Color online.) 



a ? 06? an d 0c, as shown in Fig. [TJ For completeness 
we also give the Glazer notations for our structures, 
where, e.g., a is written as a + 6°c° if consecutive oc- 
tahedra along a rotate in the same direction (in-phase) 
or a~b°c° if consecutive octahedra rotate in opposite di- 
rections (out-of-phase). 11 Multiple subscripts indicate a 
compound rotation, e.g. (j) a b denotes rotation around 
the ah axis (crystallographic direction [110] pc ). (Note 
that the 0^ are not strictly independent; however, a rig- 
orous decomposition made using irreducible representa- 
tions yields nearly identical results.) 

Biaxial strain. — The bulk crystal structure of LAO de- 
viates from the ideal Pm3m perovskite by out-of-phase 
rotations around the crystallographic [111] axis (0 a 6 C , 
Glazer system a~a~a~) that lower the symmetry to 
space group i?3c, and an accompanying rhombohedral 
distortion. We reproduce this structure by relaxing lat- 
tice parameters and atom positions in our 10-atom unit 
cell, and calculate a rhombohedral angle of 60.2° and 
an equilibrium lattice constant of 5.298 A (pseudocubic 
ao = 3.746 A), which underestimates the experimentally 
determined value by ~1%, a common artifact of LDA 
calculations. We calculate a 5 C = 5.98°, close to the ex- 
perimentally determined value. 12 

To mimic clamping to a substrate, we enforce a square 
lattice in the plane of epitaxy and constrain the value 
of in-plane lattice parameter a while relaxing the out-of- 
plane parameter c and atom positions. We also allow a 
[110] pc shear of the unit cell. The resulting structure, for 
a = 3.746 A has space group C2/c; it retains the a 5 C 
rotations of the bulk and has a relaxed monoclinic an- 
gle j3 — 90.3°. We then investigate the effect of biaxial 
strain by adjusting the value of a. For strains of -2% 
to +2% we identify two competing phases with space 
groups I A/ mem and Imma. The resulting energy ver- 
sus strain phase diagram is shown in Fig. [TJ To facil- 
itate fine-sampling near the transition we do not relax 



c in these calculations but instead maintain the equilib- 
rium volume. Later we relax this constraint and find no 
qualitative changes to the phase diagram. We find that 
the shear distortion is only energy- lowering for the C2/c 
phase and all other phases retain f3 = 90°. 

At very small strains of -0.2% to 0.1%, the ground 
state remains C2/c. As strain is increased, the C2/c 
phase becomes rapidly unstable and the system transi- 
tions to Imma or 74/mcm, depending on the sign of the 
strain. Compressive strain smaller than -0.2% stabilizes 
the 14/ mem phase, in which the lattice distortion con- 
sists of C rotations (Glazer system a°a°c _ ), whereas ten- 
sile strain greater than 0.1% stabilizes the Imma phase, 
comprised of a 5 rotations (a _ a _ c°), as shown in Fig. [3] 
Thus, the bulk- like C2/c phase is found only in a very 
narrow range of strain near the equilibrium lattice pa- 
rameter. These results are summarized in Table HI (Note 
that the precise range of stability for the three ground- 
state phases is difficult to identify from total energies. 
Thus we obtain our predicted range of the C2/c phase 
by calculating the zone-center phonons of each phase as 
a function of strain and comparing the frequencies of the 
softest non-trivial modes.) 

To examine the coupling between strain and rotations, 
we decompose the distortion from a high-symmetry par- 
ent phase with space group PA/mmm in terms of ir- 
reducible representations (or irreps) using the software 
isodisplaceJ^ We find that the dominant irreps are 
A~l and A§ , which correspond to C and a & rotations, 
respectively. The evolution of the irrep amplitudes in 
the lowest-energy structures is shown in Fig. 2J We also 
calculate C and a 5 from the angles between neighbor- 
ing octahedra to provide an approximate correspondence 
between irrep amplitude and rotation angle (Fig. [4]). 
Within I A/ mem and Imma we find a roughly linear de- 
pendence of rotation angles on lattice parameter. 

Our results indicate that the structural phases of 
strained LAO are driven by the octahedral tilting. In 
bulk LAO, compressed lattice parameters due to high 
pressure are reported to reduce rotation angles and drive 
a transition from R3c to Pm3m as a consequence of 
the greater compressibility of the A106 sites compared 
to LaOi2^~— Under biaxial strain, however, we find 
two distinct responses within the range of strains inves- 
tigated, as seen in the evolution of Al-0 bond distances 
and Al-O-Al bond angles, shown in Fig. [2j At very small 
values of strain (-0.2% to 0.1%) the changes in lattice di- 
mensions are accommodated primarily by rotations while 
Al-0 distances remain nearly constant, and the system 
remains in the C2/c phase. Under compressive (tensile) 
strain, transformation to I 4 /mem (Imma) occurs when 
0a6 (0c) goes to zero. Increased strain is then accommo- 
dated by changes to Al-0 distances, while the rotations 
remain nearly constant. 

Similar strain-induced variations in tilt patterns have 
recently been reported in a number of other perovskite 
oxides, including SrRuOsr^ BiFeOsr 19 and LaNi0 3 ^ 
The behavior of LAO is distinct from these systems, how- 
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TABLE I: Summary of the rotational modes and resulting space groups found in LAO under various strain states. 
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FIG. 3: Crystal structures of I A/ mem (left) and Imma 
(right) phases, the ground-states for compressive- and tensile- 
strained LAO, respectively. Upper figures depict [001] pc pro- 
jection (ab plane); lower figures depict [110] pc projection. La 
atoms shown in green (dark gray); Al atoms shown in blue 
(light gray) at center of oxygen octahedra. (Color online.) 



ever, in the extremely narrow window in which the parent 
phase is stable and the dramatically different structures 
stabilized by compressive and tensile strain. 

Uniaxial relaxation of strain. — We next examine the 
effect of a uniaxial relaxation of tensile strain in LAO. We 
compare the energetics of a 6 rotations (Imma phase) to 
a phase with only a rotations (space group Fmmm). 
Under biaxial strain, the Fmmm phase is ~0.5 meV 
higher in energy than Imma in the range of strains in- 
vestigated, but a uniaxial relaxation of the strain, such 
that a 7^ 6, alters the balance of energy between the two. 
Starting with atom configurations from the relaxed struc- 
tures for a = 6, we vary the ratio a/b and calculate the 
total energy for the two phases. We maintain a constant 
in-plane area ab and do not relax c. We first explore the 
system with ab = (3.746 A) 2 and find that a 0.5% distor- 
tion of a/b stabilizes a rotations relative to a & by ~1 
meV. We compare this to the system with ab = (3.85 A) 2 , 
^3% tensile strain, and find that the uniaxial relaxation 
is no longer energetically favorable and the ground-state 
retains a & rotations and a = b. 

Our results indicate that LAO films on substrates with 
a lattice mismatch of ~0.5% may, in theory, lower their 
energy through a partial relaxation such that a ^ b. 



The resulting structure exhibits a rotations (for a > b) 
in the space group Fmmm. In practice, no common 
growth substrate provides these conditions, but the re- 
sult is likely to be a general phenomenon in similar ma- 
terials. The data suggest an explanation for the experi- 
mental observation of a ^ b in some coherent films grown 
on square substrates. For example, BiFeOs films thinner 
than 50 nm grown on LAO substrates are reported to 
have a ~3.84 A and b ~3.76 A, compared to the sub- 
strate parameters a = b = 3.79 A, indicating a uniaxial 
relaxation of the epitaxial constraint^ Our results for 
LAO at 3% tensile strain, which do not predict a uni- 
axial strain relaxation, are consistent with experimental 
data for films of LAO on STO (a mismatch strain of 3%). 
Such films are reported to have square in-plane lattice 
parameters within an accuracy of 0.01 A, indicating that 
any distortion of a/b must be smaller than 0.5%*22r— 

Disabled octahedral rotations - Finally, we address the 
effect of manually disabling octahedral rotations. For this 
we treat a five-atom tetragonal unit cell in which rota- 
tion of the oxygen octahedra is forbidden by symmetry. 
At each value of epitaxial strain we relax the positions 
of the ions and calculate the electric polarization using 
the Berry phase method j 25 i 26 We find that, while the un- 
strained system is non-polar, an abrupt transition to a 
state with large polarization occurs at -3% strain. The 
polarization is out-of-plane and results from a structural 
distortion within PA/mmm symmetry in which Al and 
La displace along [001]. A compressive strain of 4% re- 
sults in a polarization of 38 jiC cm -2 relative to a cen- 
trosymmetric reference structure. These results indicate 
the existence of an incipient ferroelectric mode that is 
usually suppressed by the dominant antiferrodistortive 
rotational modes. 

Summary and discussion - In summary, we find that 
several previously unidentified phases of LAO, character- 
ized by distinct patterns of octahedral rotations, are sta- 
bilized by varying the epitaxial constraints over a range 
that is readily accessible to experiment. This remarkable 
structural softness indicates that it is unlikely that the 
bulk- type a 5 C rotations will be observed in heteroepitax- 
ial thin films of LAO on any substrate. This has impor- 
tant implications for the interpretation of structural data 
for LAO films, with particular relevance to investigations 
of the properties of LAO/STO interfaces. 
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FIG. 4: Rotational modes of LAO as a function of in-plane 
lattice parameter a: irrep amplitudes A~ and tilt angles <pj. 
(Color online.) 

We find that biaxial strain of -0.2% to 0.1% increases 
octahedral tilt angles in the bulk-like C2/c phase of 
LAO, in contrast to the experimentally observed effects of 
pressured Larger strains, however, induce a tilt-driven 
transition to one of three phases in which at least one tilt 
component is zero (74/mcm, Imma, or Fmmm). Within 
these strained phases, we find that changes in lattice pa- 
rameters related to epitaxial strain are largely accommo- 
dated by evolution of the Al-0 bond lengths, with only 
small modification of tilt angles. These findings provide 
insight to the growing list of complex oxide perovskites in 



which the rotational distortions are observed to depend 
on lattice parameters. Given the coupling between rota- 
tions and strain in LAO, and its non-polar structure, we 
propose it as a model system for exploring new non-linear 
optical techniques for probing octahedral rotations. 27 

Finally, for small values of tensile strain we find that 
a uniaxial relaxation of the strain, such that a > 6, sta- 
bilizes (j) a rotations over (j)^. Our results suggest a route 
to selectively stabilize different tilt patterns via substrate 
geometry. These results also suggest that the partial re- 
laxation of in-plane lattice parameters observed in some 
epitaxial thin films of perovskites may be driven by the 
competition between (j) a and (j)^ rotations. 
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